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Abstract 

In this we paper we recast the Cox-Ingersoll-Ross model of interest 
rates into the chaotic representation recently introduced by Hughston 
and Rafailidis. Beginning with the "squared Gaussian representation" 
of the CIR model, we find a simple expression for the fundamental 
random variable . By use of techniques from the theory of infinite 
dimensional Gaussian integration, we derive an explicit formula for 
the nth term of the Wiener chaos expansion of the CIR model, for 
n = 0,1,2, ... . We then derive a new expression for the price of a zero 
coupon bond which reveals a connection between Gaussian measures 
and Ricatti differential equations. 
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1 Introduction 



In the present paper we shall study the best known example of a term struc- 
ture with positive interest rates, namely the CIR model j2], in the context 
of the "chaotic approach" to interest rate dynamics introduced recently by 
Hughston and Rafailidis [SJ (see also £Q). By an interest rate model we mean 
the specification of a spot rate process r t and of a market price of risk process 
A t both under the "natural" or physical measure in the economy P . In the 
chaotic approach, the random nature of the model is assumed to be given by 
a probability space (Q, T ', P) equipped with a Brownian filtration (J r t)o<t<oo- 
The essence of the approach by Hughston and Rafailidis is the specification 
of the most general term structure with positive interest rates in terms of a 
single unconstrained random variable they denote X^. They then apply a 
Wiener chaos expansion to X^, and interpret the resulting terms as building 
blocks for models of increasing complexity. 

In one version of the CIR model, r t , X t are governed by the equations 

dr t = a(b - r t )dt + Cy/F t dW t , \\X t \\ 2 = X 2 r t (1) 

for some positive constants a, b, c, A with 4ab > c 2 , where Wt is a standard one 
dimensional P-Brownian motion and X t is a vector valued adapted process. 
By embedding this model inside a general class of squared Gaussian models, 
we will be led to a natural choice for the basic square integrable random 
variable associated to it. As we shall then demonstrate, the resulting 
stochastic process X t = Et[X oo] admits an explicit chaos expansion, one 
which in general includes terms of every chaotic order. 

A central idea in the present paper is the link between the Wiener chaos 
expansion and the theory of Gaussian functional integration, an essential 
tool invented to study the mathematical structure of quantum field theory. 
In fact, the mathematics underlying our example is a consequence of certain 
basic results in that theory, and can be found in for example [3]. 

The organization of the paper is as follows. In section 2 we review the 
essential ingredients for the construction of positive interest rates models in 
both Flesaker-Hughston and the state price density approaches, and then 
compare these approaches to the recently introduced chaotic representation 
of Hughston and Rafailidis [3] . We end the section by describing the structure 
of the Wiener-Ito chaos expansion and show how it can be expressed in terms 
of a certain generating functional acting on the space L 2 (M. N ). 
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In section 3, we describe the squared Gaussian formulation of the CIR 
model and show the spot rate process can be explicitly computed. Based on 
this representation, we state the form of the random variable X^, and give 
a proof that it lies in L 2 . In section 4, we state the exponential quadratic 
formula which is the main technical tool in this paper. It is a formula for the 
generating functional of random variables of the form X = e~ Y for Y lying 
in a general class of elements in the second chaos space 7^2- In section 5, 
we compute the generating function for the random variables X t in the CIR 
model and, as the main result of the paper, derive their chaos expansion. In 
section 6, we show that the usual CIR bond pricing formula has a natural 
derivation within the chaotic framework. 

Three appendices focus on the theory of Gaussian functional integration 
and its relation to the Wiener chaos expansion. Appendix A explores the 
white noise calculus. Appendix B states and provides a proof of the gener- 
ating functional theorem. Appendix C provides a proof of the exponential 
quadratic formula. 



2 Positive Interest Rates 



2.1 State price density and the potential approach 

Rather than focus on the spot rate process, one can model the system of 
bond prices directly. Let P t T, < t < T denote the price at time t for a 
zero coupon bond which pays one unit of currency at its maturity T. Clearly 
Pu = 1 for all < t < oo and furthermore, positivity of the interest rate is 
equivalent to having 

Pts<Ptu, (2) 

for all < t < u < s. 

A general way to model bond prices [121 E] is to write 

P. = 3™ (3) 

for a positive adapted continuous process V t , called the state price density. 
Positivity of the interest rates is then equivalent to V t being a supermartin- 
gale. In order to match the initial term structure, this supermartingale needs 
to be chosen so that i?[Vr] = -Pot- If we further impose that P t -> as 
T — > oo, then Vt satisfies all the properties of what is known in probability 
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theory as a potential (namely, a positive supermartingale with expected value 
going to zero at infinity). 

It follows from the Doob-Meyer decomposition that any continuous po- 
tential satisfying 



E sup V t 2 < oo 



0<t<oo 



(4) 



can be written as 

Vt = E t [A 00 ]-A t , (5) 

for a unique (up to indistinguishability) adapted continuous increasing pro- 
cess A t with E(A^ I0 ) < oo. Therefore, the model is completely specified by 
the process A t , which can be freely chosen apart from the constraint that 



E 



~dAr 



' dT ' 



(6) 



2.2 Related quantities and absence of arbitrage 

An earlier framework for positive interest rates was introduced by Flesaker 
and Hughston [3J eq. (8)], who observed that any arbitrage free system of 
zero coupon bond prices has the form 

C° h s M ts ds , . 

Pit = L for < t < T < oo. (7) 

J t h s M ts ds 

Here hx = is a positive deterministic function obtained from the initial 

term structure and M ts is a family of strictly positive continuous martingales 
satisfying M 0s = 1. Any such system of prices can be put into a potential 
form by setting 

/oo 
h s M ts ds. (8) 

The converse result is less direct and was first established by Jin and Glasser- 
man lemma 1]. 

These equivalent ways of modelling positive interest rates can now be 
related to other standard financial objects. A particularly straightforward 
path is to follow [121 proposition 1] : given a strictly positive supermartingale 
V t , there exists a unique strictly positive (local) martingale A t such that the 
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process B t = A t /V t is strictly increasing. We identify B t with a riskless 
money market account initialized at B — 1 and write it as 

B t = exp( / r s ds), (9) 
Jo 

for an adapted process r s > 0, the short rate process. 

A sufficient condition for an arbitrage free bond price structure in the 
potential approach is to require that the local martingale A t be in fact a 
martingale, since it can then be used as the density for an equivalent mar- 
tingale measure. It is an interesting open question in the theory to isolate 
what conditions on the potential V t would suffice for that. 

The formulation up to this point is quite general, in the sense that is does 
not make use of any particular structure of the underlying filtration (other 
than the usual conditions). Let us now assume that Tt is actually generated 
by an A-dimensional Brownian motion Wt- The market price of risk then 
arises as the adapted vector valued process X t such that 

dA t = -\ t A t dW t . (10) 

It is also immediate to see that the state price density process is the solution 
to 

dV t = -r t V t dt - X t V t dW t , (11) 

so that the specification of the process V t is enough to produce both the short 
rate and the market price of risk At. 

It had already been remarked by Flesaker and Hughston j^j that in the 
Brownian filtration with finite time horizon any positive interest rate model 
in their formulation corresponds to a model in the HJM family with positive 
instantaneous forward rates ftr- The converse result that any interest rate 
model in HJM form with positive instantaneous forward rates can be written 
in the Flesaker-Hughston form was also obtained by Jin and Glasserman 
theorem 5]. In order to prove this result they found a rather technical 
necessary and sufficient condition for positivity in terms of the volatility 
structure of the HJM form, confirming that the HJM formulation is not the 
most natural one to investigate positive interest rates. 

2.3 The Chaotic Approach 

We have seen in the potential approach that the fundamental ingredient to 
model the random behaviour of the interest rates is the increasing process A t 
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in the decomposition V t = E^Aoo] — A t , whereas in the Flesaker-Hughston 
construction the corresponding role is played by the martingales M ts . 

In jS], Hughston and Rafailidis introduced an elegant construction of 
general positive interest rate models based on a Brownian filtration using 
simpler fundamentals. Assuming that the state price density V t is a potential 
satisfying 



E 



r x V,ds 



< oo, 



(12) 



then integrating (jllj) on the interval (t,T), taking conditional expectations 
at time t and the limit T — > oo, one finds that 



V t = E t 



r,V,ds 



Now let a t be a vector valued process such that 

Of (J t = r t V t , 
and define the square integrable random variable 



(J,dW,. 



It then follows from the Ito isometry that 

V^E^X'J-E^] 2 , 



(13) 



(14) 



(15) 



(16) 



which is called the conditional variance representation of the state price 
density Vt- To obtain the connection between this representation and the 
Flesaker-Hughston framework, observe that a direct comparison between 
(fHl) and © gives that 



Is 



E f 



Uc • a. 



(17) 



Similarly, by comparing the conditional variance representation (|16j) with the 
decomposition (jSJ), we see that 



where X t = E^X^ 
decomposition that 



EtiX^-Xf^EtiA^-At, 

It follows from the uniqueness of the Doob-Meyer 
A t = [X,X] t , 
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that is, the quadratic variation of the process X t . 

Conversely, given a zero- mean random variable X^ G L 2 (Q, J 7 , P), the 
representation (fTBj) defines a potential V t , which can then be used as a state 
price density to obtain a system of bond prices. The issue of absence of arbi- 
trage can then be addressed in terms of necessary and sufficient conditions on 
Xqo, and is by and large an open question at this point. The construction in 
[S] flows in the opposite direction, in the sense that the authors first enumer- 
ate a series of axioms to be satisfied by an arbitrage-free interest rate model 
and then obtain a square integrable random variable I m corresponding to 
it. 

2.4 Wiener chaos 

As Hughston and Rafailidis also observed, the L 2 condition on X^ is nec- 
essary and sufficient for I ro to have the type of orthogonal decomposition 
known as a Wiener chaos expansion j^j. They interpret the different or- 
ders of this decomposition as basic building blocks for models of increasing 
complexity. 

Let W t be an iV-dimensional Brownian motion on the filtered probability 
space (ft, J 7 , {jF t } tgR+ , P). We introduce a compact notation 



s,n) G A = R + x {!,..., 



N} 



and express integrals as 




(18) 



For each n > 0, let 



H n {x) = (-ire* 2 / 2 ^e- 2 / 2 (19) 

be the nth Hermite polynomial. For h G L 2 (A), let \\h\\ 2 = j^h(r) 2 dT and 
let W(h) be the Gaussian random variable J A h{r)dW T . The spaces 

H n = spm{H n (W(h))\h G L 2 (A)}, n > 1, 
Ho = C 
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form an orthogonal decomposition of the space L 2 (Q, Too, P) of square inte- 
grable random variables: 

L 2 (Q,T oo ,P) = (B^ =0 n n 
Each H n can be understood completely via the isometries 

J n : L 2 (A n ) — > Ti n 

given by 

/„-./„(/„)=/ Ufa , . . . , T n )dW Tl . . . dW Tn (20) 

where A n = {fa, . . . ,T n )|r, = (s^/i,) G A, < s 1 < s 2 < ■ ■ ■ < s n < oo}. 

With these ingredients, one is then led to the result that any X G 
L 2 (Q, Too, P) can be represented as a Wiener chaos expansion 



X = Y J UU) 



(21) 



n=0 



where the deterministic functions /„ G L 2 (A n ) are uniquely determined by 
the random variable X. 

A special example arises by noting that for h G L 2 (A) 



n\J n (h® n ) 



nrr ( W{h) \ 

n V \\H ) 



(22) 



and furthermore 
exp 



In the notation of quantum field theory (see Appendix A), this example 
defines the Wick ordered exponential and Wick powers 



: exp[W(h)\ : 



exp 



W(h) 



h(r) 2 dr 



:W(h) n : = n\J n (h® n ) 



(24) 



Generating functionals provide one systematic approach to developing 
explicit formulas for the terms of the chaos expansion in specific examples. 
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Theorem 2.1 For any random variable X G L 2 (fl, J 7 ^, P) , the generating 
functional Z x (h) : L 2 (A) — > C defined by 



Zx{h) = E 



X exp 



W{h) -\f h(r) 2 dr 



(25) 



is an entire analytic functional of h G L 2 (A) and hence has an absolutely 
convergent expansion 

Z x (h) = J2F { x ) (h) (26) 



n>0 



where 



F [ x\h) 



fx\n, T n )h(Ti) . . . h(r n )dTi . . . dr n . 



(27) 



The n-th Frechet derivative of Z x at h = 0, fx\ T i> ■ ■ • > T n), lies in L 2 (A). 
Finally, the Wiener-Ito chaos expansion of X is 



n>0 

Proof: See appendix B. 



X = J2[ f i x\r 1 ,...,r n )dW Tl ...dW Tt 



(28) 



□ 



3 Squared Gaussian models 

A number of authors |H],|S], [H] have observed that the CIR model [2j with an 
integer constraint N = ^ G N + \ {0, 1} lies in the class of so-called squared 
Gaussian models. By introducing an IR^-valued Ornstein-Uhlenbeck process 
R t , governed by the stochastic differential equation 

dRt = -^Rt dt + ^dWt (29) 

where W t is iV-dimensional Brownian motion, the Ito formula together with 
Levy's criterion for Brownian motion shows that the square r t = R} t Rt satis- 
fies © where W t = jl{R\R t )- 1/2 RfdW t is itself a one-dimensional Brownian 
motion. Note that here, and for the remainder of the paper, we suppress vec- 
tor indices by adopting a matrix multiplication convention, including f for 
transpose, in which for example R\^{t)dW t denotes Ylu, v =i Rtl^ifydW". 
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We focus on a general family of interest models which includes this exam- 
ple, the so-called extended CIR model and more. Note that we work always 
in the physical measure and thus to specify the term structure model one 
needs to determine the market price of risk vector X t as well as the spot rate 
process r t . 

Definition 3.1 A pair (r t ,\t) of (O, J 7 , Tt, P) processes is called an N- 
dimensional squared Gaussian model of interest rates (N >2) if there is an 
M iV -valued Ornstein-Uhlenbeck process such that r t = R\Rt and X t = XR t . 
R t satisfies 



where a, 7, A are symmetric matrix valued and R vector valued determinis- 
tic measurable functions on M. + . W is standard N -dimensional Brownian 
motion. In addition we impose boundedness conditions: 

• there is some constant M > such that a(t) > M and |A(t)|~ 2 > M 
for all t 

The exact solution of (j5U|) is easily seen to be 



dR 



\ = a(t)(R(t) - R t )dt + y(t)dW t , R 



t=o — 



Rq 



(30) 




(31) 



where 




(32) 



and K(t, s), t > s is the matrix valued solution of 




a(t)K(t,s) 0<s<t 
< t 



(33) 



which generates the Ornstein-Uhlenbeck semigroup. 
By pip, the state price density process is 




(34) 



We thus have a natural candidate for the random variable X, 



00 ■ 




(35) 
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where the E -valued process 



a t = exp 



R\ ( - + — ) R s ds + -Rt\dW s 



Rt. 



(36) 



is the natural solution of a]a t = r t V t . 

Before proceeding to analyse X^ in detail we show that X^ is square 
integrable. 

Proposition 3.2 Iwht^,^ E[X^\ = 1 

Proof: By the Ito isometry and Fubini's theorem for the Ito integral, 

f-T 



E[X, 



E 



at atdt 



o 



(37) 



From a\a t dt = r t V t dt = -dV t - V t \UW t it follows that E[X%] =E[1-V T }. 
We observe that this is true for any pair X t , V t defined using the formal 
relations (fTT j) .(fH |) .(fT5 J) of the chaotic approach. All we need to prove now is 
that the supermartingale V t is indeed a potential when (r t , X t ) is a squared 
Gaussian model. 

For any < e < 2M, V T = e~ Y ^ where Y l = ^ R\{l - ^)R t dt is 
positive and 

Y 2 = -i- ( T \\r\{\ + efX 2 R t dt + (1 + e)RlXUW t 

By the Holder inequality 

E[V T ] < {E\e^ +1 l^]) el{1+e) {E[e^ 1+ ^}) 1/(1+e) 

with the second factor equal to 1, since e~( 1+e ) y2 is an exponential martingale. 
Now Yi is a positive random variable for which a direct computation shows 

mean(Yi) = CiNT (1 + 0{T)) (38) 
var(Y 1 ) = C 2 NT(l + 0{T)) (39) 

for positive constants Ci, C 2 . An easy application of Chebyshev's inequality 



Prob Yi < 



dNT 



< o{ 



NT' 



(40) 



then implies that lim^oo -E^Vr] = 0. 



□ 
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4 Exponentiated second chaos 



The chaos expansion we seek for the CIR model will be derived from a closed 
formula for expectations of e~ Y for elements 

Y = A+ [ B( Tl )dW Tl + [ C{r U T 2 )dW T1 dW T2 (41) 

J A J A 2 

in a certain subset C + C TL<2 = 71q ®71\@ 71.2 ■ in the integrals above, 
recall that compact notation using r's carries a summation over vector in- 
dices as well as integration over time. The formula we present is well known 
in the theory of Gaussian functional integration jH Chapter 9]. In prob- 
ability theory, this result gives the Laplace transform of a general class of 
quadratic functionals of Brownian motion. Many special cases of this result 
have been studied in probability theory, see for example [TSJ Chapter 2] and 
the references contained therein. 

If in (}4~Tj) we define C(ti,t 2 ) = C(t 2 ,ti) when t\ > r 2 , then C is the 
kernel of a symmetric integral operator on L 2 (A): 

POO 

[Cf\(r)= / Cfard/CrOdTi (42) 
Jo 

Recall that Hilbert-Schmidt operators on L 2 (A) are finite norm operators 
under the norm: 



\\C\\ 2 HS = l C(T l ,T 2 ) 2 dT l dT2 
J A 2 



We say that Y G 71<2 is in C + if C is the kernel of a symmetric Hilbert- 
Schmidt operator on L 2 (A) such that (1 + C) has positive spectrum. 

Proposition 4.1 Let Y e C + . Then 
E[e~ Y ] = [det 2 (l + C)n 1/2 



exp 



-A + - 
2 



If S(r 1 )(l + C7)- 1 (r 1 ,r 2 )S(r 2 )dr 1 dr 2 
z JAo 



(43) 



Remark 4.2 The Carleman-Fredholm determinant is defined as the exten- 
sion of the formula 

det 2 (l + C) = det(l + C) exp[-Tr(C)] (44) 

from finite rank operators to bounded Hilbert-Schmidt operators; the operator 
kernel (1 + C) _1 (ri, r 2 ) is also the natural extension from the finite rank case. 
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Proof: See Appendix C. 



□ 



Using this proposition, it is possible to deduce the chaos expansion of the 
random variables e~ Y , Y G C + , a result known in quantum field theory as 
Wick's Theorem: 

Corollary 4.3 IfY = J A ^C(ri,r 2 )dW Tl dW T2 E C + , then the random vari- 
able X = e~ Y has Wiener chaos coefficient functions 



/ n (ri, . . . ,r n ) 







n even 
n odd 



where K = [det2(l + C)j and for n even, Q n is the set o/Feynman graphs 
on the n marked points {n, . . . ,r n }. Each Feynman graph G is a disjoint 
union of unordered pairs g = (r gi , r 92 ) with U 9e c g = {ri, . . . , r n }. 

Proof: The generating functional for X = e~ Y is 



Z x {h) 



E 
E 



Xexp h(r)dW T -~J h(r) 2 dr 
exp (^J h(r)dW T -^J h(r) 2 dr - J C(t u T 2 )dW T1 dW T . 



so we can use Proposition 14.11 with A = | J h(r) 2 dr and B(r) = —h(r) 
which yields 



Z x (h) 



det 2 (l + C)- 1/2 
1 



exp 



hKnMn,^) - (1 + C)- 1 (T U T 2 )}h(r 2 )dT 1 dr 2 



A 2 



.(45) 



Using the last part of Theorem 12. 1| the result comes by evaluating the nth 
Frechet derivative at h — 0, or equivalently by expanding the exponential 
and symmetrizing over the points Ti, . . . ,T n in the n/2th term. 

□ 
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5 The chaotic expansion for squared Gaus- 
sian models 

We now derive the chaos expansion for the squared Gaussian model defined 
by (|30j). In view of (|35|) it will be enough to find the chaos expansion for cr^, 
T < oo. We start by finding its generating functional Z a ^. For h,k G L 2 (A), 
define the auxiliary functional Z(h, k) = E [e _yT ] with 



Yt 



R 



^2 

, i - H ) Rtdt . 

1 2 4 / 2 



T 



RjXdWt 



h\t)dW t 



1 

+ 2 



Jo 

T i-T 

h\t)h(t)dt- / k\t)R t dt. 
n Jo 



(46) 



Proposition 5.1 Z(/i, /c) an entire analytic functional on L 2 (A) x L 2 (A). 
Moreover 

5Z(/i, k) 



lim 



fc=0 



(47) 



where Z a ^{h) is defined by 1^3)) twi/i X = cr^,/i 



iV. 



Proof: Analyticity in (h, k) follows by repeating the argument given in Ap- 
pendix B. By the definition of Frechet differentiation and continuity of the 
t -> T~ limit, g7j) follows. 

□ 

We want to use Proposition ^. II in order to compute Z(h, k). Substitution 
of (jSH), into the first term of fljB)) leads to 



1 A 2 
2- + -4 



R\ ( ~ + '— ) fltd* 



r 



7 (t)dWi 



+ / 7(*i 

'A 2 

+ / tr< > 7(t) 



£ K T (ti,s) (l + ^pjK T (s,t 2 )ds 7 (t 2 )dW tl dW t 
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where we define ifr^ij^) = < T)K(ti,t 2 ). For the second and the last 
terms of (1461) we have 



R\\dW t 



2 Jo 



R\t)\dW t + - I tr 



2 Jo 



^{s)K T (s,t)Xds 



dt 



+ 



tf(t)Rtdt 



o 



7(*i)l4(*i> *2)A + AK T (ti, t 2 )i(t 2 ) ) dW tl dW t2 , 

T / pT 



tf{t)R(t)dt 



o 



tf(s)K T ( Sj t)ds 7(t)dW( 



'0 wo 

Thus the exponent Yr appearing in (|4fi|) has the form of (|4~T|) with 



An 



j T (t) £ KUt,s)(~ + ^pjK T {s,t)ds 7 (t)Jdt 



+ / tr 7 (t) 



Br(t) 



-h{t)-j{t) [ Kl(t,s)k(s)ds + l\R(t) 
Jo 2 



T 



+ 7 (t) / + 



A s 



C T (ti,t 2 ) 



7(*i 
1 



7 (ti)i4(*i. *2)A + A\RT r (ti, f 2)7(^2 



7& 



It is clear that the operator CV has Hilbert-Schmidt norm ||Ct||# 5 



(9(T). Moreover, if we denote by ^K T (1 + — J iiVy, 7-K^A and XKt'J the 
operators whose kernels appear in the expression above, then Ct can be 
written as 



C T = 7 ^t7 + ^(7^A + 1)(A^t7+1)-^ 



(48) 



from which we see that (1+CV) is positive. Therefore, we can use Proposition 
14.11 for E[e~ YT ], leading to a general formula for the generating functional 
Z(h,k): 
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Z(h, k) = det 2 (l + CV)~ 1/2 exp (-~trC T ^J 



x exp 
x exp 



dt 



X 



A 2 

2 

T 



K T {sM)l{ti)ds 



(1 + C T )- L (t u t 2 ) 



h(t 2 ) + 7 (* 2 ) / K J T (t 2 , s)k s ds - -XR(t 2 ) 



-l(t 2 )J K^ T (t 2 ,s)(l + ?p}R{s)d S 
Differentiation once with respect to k then yields 

cT V 



dtidtr 



(49) 



Z aT {ti) = M T exp 







^ t(t) G + t) A (*)+^ f (W*) 



(ft 



x <^ - J R + J fr T7 (l + Cr)- 1 



Ai? 

T 



/i - ^ - 7_?4 ( 1 + — ) i? 



A 2 
2 



x exp <; - 



A 2 



/it 



x 



Ai? 



(t 2 )dtidt 2 



A 2 



(T) 

(Oa + Cr)- 1 ^,^) 
(50) 



where 



M T = e-^ trCT (det 2 (l + C T ))" 1/2 = (det(l + C T )y 1/2 . (51) 

By comparing (}4*9*|) and ()50|l. the reader can observe our use of an operator 
notation which suppresses some time integrals, for example in the very last 
term 



7 l4 1 1 + '— ) R 



A 2 
2 



(t)= 7 (*) / KUt,s)[l + ^)R(s)ds, 



and similarly for other terms. 
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These formulas simplify considerably if the function R vanishes, which 
is true in the simple CIR model of ()29|) when r = 0. In this case we have 
a{t) = a/2 and j(t) — c/2 (here certain scalars are to be understood as 
multiples of the identity matrix), so that K T (s,t) = e~ a ( s ~ t " 2 l(t < s < T) 
and 

Crit.M) = ^(i + f) [e-fl tl - t2 l-ei^+ t2 - 2T )]+|Ae-fl* 1 - fe l(52) 

Moreover, the previous expression for Z aT (h) reduces to 
Z aT (h) = M T [K r7 (l + C T )- 1 /i](T) 



cxp 



-~ [ rf(t)h(t)dt+~ f h\h)(l + C T )-\t^t 2 )Kt2)dt x dh 
1 Jo 1 Ja 2 



We can then easily evaluate the nth Frechet derivative of Z UT at h — 
as in the proof of Corollary 14.31 and determine the following partly explicit 
form for the nth term of the chaos expansion. 

Theorem 5.2 The nth term of the chaos expansion of a? for the CIR model 
with initial condition r = is zero for n even. For n odd, the kernel of the 
expansion is the function f^\-) : A n — > W 

f T (t 1 ,...,t n )=M T U L ^) (53) 
Geg* geG 

where 

( [cva+cvr 1 ]^,^) Tig 

L{9) = (54) 
[ {K Tl {l + C T y l ){T,t g2 ) Teg 

Here, Q* is the set of Feynman graphs, each Feynman graph G being a par- 
tition of {t\, . . . ,t n ,T} into pairs g = (t 9l ,t g2 ). 

The chaos expansion for Xoo itself is exactly the same, except that the 
variable T is treated as an additional Ito integration variable. The explicit 
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expansion up to fourth order is: 
A\ = / M T [K T7 (l + C T )- 1 ](T,t 1 )« 1 rfW T 



+ / M T [K T7 (l + C , T )- 1 ](T,t 3 )[C , T (l + C' T )- 1 ](t 1 ,t 2 )^ 1 « 2 ^3^ T 
Ja 4 

+ / M T [K Tl (l + C T )- 1 ](T,t 2 )[C T (l + C T y 1 }(t 1 ,t 3 )dW tl dW t2 dW t3 dW T 
M T [K Tl (l + C f T )- 1 ](T,t 1 )[C T (l + C T )~ 1 ](t 2 ,t 3 )dW tl dW t2 dW t3 dW T 



(55) 



6 Bond pricing formula 



In this section we give a derivation of the price of a zero coupon bond in the 
CIR model. Recall from section |2] that these are given by 



P tT = E t [V t ~ l V T ]. 



(56) 



To keep things "as simple as possible, but not any simpler", we take A = 
so Vt = exp[— J Q r s ds], or in terms of the squared Gaussian formulation, 



(57) 



Vt = exp — / w s R s ds 

L Jo 

As we have seen in the previous section, for t < s <T 

R» = K T (s, t)R? + ^J' K T {s, Sl )dW£ 

rp 

hence — log[Vp Vr] = Ji (Rs) 2 ds can be written as 

-T r s 2 



(58) 



Jt 



4 4R f r r S2 

-(R?) 2 c T (t,t) + / c T (t,s)dw s ^ + 2 / c T ( Sl , S2 )dw s yw s ^ 

„ l c c Jt Jt Jt 

+N J C T (s,s)ds 



(59) 



where Ct(s\, s 2 ) is given by (J52|) with A = 0. 



Taking the conditional expectation of V t 1 Vt by use of Proposition 14.11 
leads to the desired formula 



P tT = [det(l + 2C T )] 



-N/2 



Hexp --(Rtf(C T (l + 2C T )- 1 )(t,t) 



(60) 
(61) 



Thus PtT has the exponential affine form exp[— (3(t, T)r t — a(t, T)] with 

P(t,T) = ±[C T (l + 2C T y 1 }(t } t) 
c z 

N 



-log[det(l + 2C r )] 



a(t,T) 

The known formula has the same form, with 

2(^ r -*) - 1) 



(62) 



P(t,T) 
a{t,T) 



(p + a)(e^*) - 1) + 2p 5 



a 2 + 2c 2 



2ab, 
— log 

P 2 



2pe 



(a+p)(T-t)/2 



_(p + a)(e^ T -*) - 1) + 2p_ 



(63) 



which can be derived as solutions of the pair of Ricatti ordinary differential 
equations 



dt 
da 

~dt 



c^5 2 
2 

-abp 



+ a(3-l 



(64) 



One can demonstrate using power series expansions that (jnzjl do in fact 
solve the Ricatti equations and hence agree with the usual formula. This 
example points to the rather subtle general relationship between kernels such 
as (1 + 2Ct)~ 1 and solutions of Ricatti equations deserving of further study. 



7 Discussion 

We have shown how the CIR model, at least in integer dimensions, can be 
viewed within the chaos framework of Hughston and Rafailidis as arising 
from a somewhat special random variable X^. This random variable can 
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be understood as derived from exponentiated second chaos random variables 
e~ Y ,Y G C + ). Such exponentiated C + variables form a rich and natural 
family which is likely to include many more candidates for applicable interest 
rate models. Although their analytic properties are complicated, there do 
exist approximation schemes which can in principle be the basis for numerical 
methods. 

On the theoretical side, this family is distinguished by its natural invari- 
ance properties. Most notably, as will be investigated elsewhere, it is invari- 
ant under conditional JF 4 -expectations: logE t [e~ Y ] G C + whenever Y G C + . 
Note in particular that this implies that these can be used as the Radon- 
Nikodym derivatives of measure changes which generalize the Girsanov trans- 
form, and which can greatly enrich the tools applicable in finance. 

This application of the chaos expansion to squared Gaussian models also 
illustrates a deep connection between methods developed for quantum field 
theory and the methods of Malliavin calculus. Many of the very rich ana- 
lytic properties of this example reflect well known techniques widely used in 
mathematical physics. 

To conclude, even if the representation of the CIR model and its gen- 
eralizations we present is not simple, it does show the way to the use of a 
powerful set of analytical techniques which might prove to be very useful in 
financial modelling. 

A White noise calculus 

Here we describe the white noise calculus, which can be regarded as a refor- 
mulation of the calculus of Wiener measure into concepts familiar to practi- 
tioners in quantum field theory as Gaussian functional integration. We follow 
the discussion of 0ksendal [TUj . 

Let S be the Schwartz space of smooth functions on M + of rapid decrease, 
and S' its topological dual, the space of tempered distributions on R + . If 
(j) G S' and / G S, we write 



for their canonical pairing. We also use the formal notation 
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to represent the 'smearing' of the distribution over the test function /. It 
acquires a rigorous meaning, however, in the cases where is itself a function 
on R + for which the pointwise product </>(s)f(s) is integrable for all / e S. 
Define now the functional 

S{f} = e -5ll/ll 2 =e ~\(fJ) L ^ ( 65 ) 

where (•, -) L 2 denotes the real-valued inner product in L 2 (IR + ). Observe that 
this functional satisfies 

1. (continuity) /„ — > / in S implies that S{f n } — > S{f} ; 

2. (positive definiteness) Yldj=i^i c j^{fi ~ fj} ^ f° r an fi^S an d all 
Ci e C ; 

3. (normalization) S^O} = 1. 

It follows from the Bochner-Minlos theorem that there exists a unique 
Borel probability measure \i on S' such that S{f} corresponds to a moment 
generating functional, that is, for all / e 5 

/ eWf)dn(4>) = S{f} = e-ll^l 2 / 2 . (66) 

The measure space (S', £>, //), where i3 is the Borel a-algebra of S', is called 
the white noise probability measure. 

In Euclidean quantum field theory, a given Borel measure \x on S' charac- 
terizes the family of 'fields' G S' through the properties of the random vari- 
ables </>(f) : S' — > R obtained for each f E S. One tries to construct measures 
pi so that the generating functional S{f} satisfy the so-called Osterwalder- 
Schrader axioms, in order to guarantee that the fields have certain required 
physical properties. The family of Euclidean free fields is obtained when 

Sc{f} = e-^= [ e^// c (0), 
Js' 

where C is the integral kernel of a positive, continuous, nondegenerate Eu- 
clidean covariant bilinear form C on S x S. We see that the special case 
S{f} = e-H-fll 2 / 2 is obtained when C(s,t) = 5(t - s), called the "ultralo- 
cal" covariance. For each e S', the random variables {4>{f) : / 6 5} C 
L 2 (S', B, n)} form a Gaussian family with mean zero and covariances 

EMf)*(9)\ = I f(s)g(t)C(s,t)dsdt=(f,g). (67) 

JR + xR + 
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The theory of martingales makes its appearance in white noise calculus 
through the concept of Wick ordered random variables. To introduce them, 
recall that the Hermite polynomials are defined by 

H n (x) = (-lTe*£;(e-£), 
from which it is easy to see that they satisfy the recurrence relation 



X ~ dx ' Hn -^ X ' 



= H n (x), n = l,2,... 
H (x) = 1. 

The first couple of Hermite polynomials are then 

H (x) = 1 

Hi(x) = x 

H 2 (x) = x 2 -l 

H 3 (x) = x 3 -3x,... 

We begin by defining the Wick ordered exponential for any / e S to be 

. e m := e <t>u)-\\f\\v\ (68) 

Then we have 



Proposition A.l For any f e S, 



: e T 



-H r 



n>0 



where H n is the nth Hermite polynomial. Moreover, for any f,g e S, 

(f,g) 



\m) Hm vw) 



\9\ 



(69) 



(70) 



Proof: For any a, b G R we have the absolutely convergent expansions 

y-a/b) 2 /2+a 2 /(2b 2 ) 



e a~b 2 /2 = e -( 6 -a/ 6 )2/2+a 2 /(26 2 ) 



V— f-ire a2/(2b2) — (e~ x2/2 ) 
^ nV ' dx n \ J 



n>0 



x=a/b 



n>0 U - 



(71) 
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where the last line makes use of the defining property of the Hermite poly- 
nomials. Using this with a = <j)(f),b = \\f\\ gives (JoTJj) . The orthogonality 
relation ()70|) follows by expanding the identity 

E^. e^ f) :: e<^ } :] = e {f ' 9) (72) 

in powers of /, g and comparing to the expansion derived from (JoT?|) • 

□ 



From this we define the Wick ordered monomials as the random variables 

(73) 



so that linearity and convergent power series imply that 

. e m . y ■ »(/)" ■ 

n>0 

Formally, we express the Wick ordered monomials as 

■ 4>(fT ■= / f(si) ■ ■ .f(s n ) : (p si . . . 4> Sn : ds x ...ds n 
and from the orthogonalization (fTUj) we can define the Wick products 



:o(f.)...o(f„) 



/ . . . fn(s n ) ■ 0si • • ■ 4>s n ■ dsx... ds n ,(7A) 



7i n is defined to be the span of {: <f)(f) n '■ \f G <S} and consists of pre- 
cisely the random variables (n!)" 1 J Rn /(si, . . . , s n ) : <f) si . . . <f) Sn : ds\ . . . ds n 

where / lies in T n the L 2 completion of the space of symmetric functions in 
iS® n . These subspaces form an orthogonal decomposition of L 2 (S',B,fi). In 
quantum field theory this is known as the Fock space decomposition of the 
Hilbert space of quantum states into n particle sectors for n > . 

The relationship between Wiener measure and the white noise measure 
is to identify </>(f) = J /(s)0 s rfs with W(f) = J f(s)dW s for all / G S C 
L 2 (IR + ). One is then lead to the formal relation S = dW s /ds, that is "white 
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noise" is the "derivative" of Brownian motion. This identification extends to 
all orders in the chaos expansion via: 

(niy 1 / f(si, . . . , s n ) : 4> S1 . . . 4> Sn : ds x ... ds n = / /(si, . . . , s n )dW Sl . . . dW s . 

(75) 

where the bijection /<->•/ between T n and L 2 (A n ) is the restriction map and 
its inverse. Finally, this leads to the identification L 2 (S', B, fi) = L 2 (Q, JF^, P). 

One useful consequence of Wick products and the chaos expansion is that 
conditional expectations can be handled systematically. 

Proposition A. 2 Let f e L 2 (R + ) and t G [0, oo). Then 

E[. e m : \F t ] = : e^)) : 
E[: (cf>(f)r : \F t ) = :(0(/ W ))": (76) 

where [/(t)](s) = l(s < t)f(s). In other words, the processes : e^ (t) ) :, 
: {4>{f(t))) n '■ are martingales. 



B Generating functional for chaos coefficients 

We derive Theorem 12 .11 in the notation of the white noise calculus. To begin, 
we recall the definition of analyticity for a function between complex Banach 
spaces. 

Definition B.l Let f : B — > C where B,C are complex Banach spaces. 
Then f is analytic at a point x £ B if 

1. for any finite set {x±, . . . , x n } C B of points the function F : C n — ► C 

F(Ci,...,C»)=/(z + X;CiSi) (77) 

i 

is analytic at £ C n ; 

2. f is continuous at x. 

Proof of Theorem 2.1 One can check directly that the map h i— >: e^ h ^ : is 
entire analytic between L^(M + ) and L^S', £>, fi). Similarly, Z x {h) is analytic 
for any h £ L^(R + ). 
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Finally, to verify (J28J), it is enough to take the expectation of the equation 
multiplied by : <f)(g) n : for arbitrary g G L 2 (A), n > 0: 

cT 



E[X : ct>{g) n '] 



whereas by (|70*j). 



-^x(A^) 



A=0 



™ ! / /x ) ( s b---) s n)5 , ( s i)---fi , ( s n)dsi...(is n (78) 



^[^(m!)- 1 / /^ m) ( Sl , ...,s m ): <t>( Sl ) . . . (j>{s m ) :ds 1 ... ds m : <f>(g) n :] 
m >o 



n! / g(s 1 )...g(s n )f i x\s 1 ,...,s n )ds 1 ...ds 



(79) 
□ 



C The exponential quadratic formula 

Here we prove the formula (|43jl in the context of one-dimensional white noise 
calculus. The proof extends easily to the multidimensional case. 

Proof of Proposition I4.lt 

1. (Step 1) First we note that Y e C + can be approximated in the Hilbert 
space 7i<2 by random variables of the form 

M 

Y = A + J2[bi<f>{9i) + <k : 4>( 9l ) 2 : /2] (80) 

i=l 

with {(yfj} a finite orthonormal set in L 2 (A) and numbers q > — 
Then Xj = 0(#j) form a collection of independent iV(0, 1) random 
variables. For K of this type, there is a factorization 

E[e- Y ] = e- A {27i)- M ' 2 \\ [ exp[-6 4 x - (1 + a)x 2 /2 + Cl /2)]dx (81) 

i 

into one dimensional Gaussian integrals. Each integral gives the factor 
(2tt) 1 / 2 (1 + c,)- 1/2 exp i[c, + 6 2 (1 + q)" 1 ] (82) 
25 



leading to a formula for E[e Y ] which agrees with (j4.1j) for Y of this 
form. 

2. (Step 2) Since the formula is true for Y in a dense subset of C + , it is 
now enough to prove that the map Y i— > i?[e _y ] is continuous in C 
provided the kernel C satisfies the stated conditions. By the definition 
of the Carleman-Fredholm determinant 

det 2 (l + C)= exp[Tr[log(l + C) - C). (83) 

Since log(l + x) — x = 0(x 2 ) for x — > we see that det 2 (l + C) is well- 
defined and continuous for C > — 1 and Hilbert-Schmidt. Therefore, 
the entire right hand side of (|4.1|) is continuous in C for Y G C + . 

□ 
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